Magnetic and Electronic Properties of Metal- Atom Adsorbed Graphene 
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We systematically investigate the magnetic and electronic properties of graphene adsorbed with 
diluted 3d-transition and noble metal atoms using first principles calculation methods. We find 
that most transition metal atoms (i.e. Sc, Ti, V, Mn, Fe) favor the hollow adsorption site, and the 
interaction between magnetic adatoms and vr-orbital of graphene induces sizable exchange field and 
Rashba spin-orbit coupling, which together open a nontrivial bulk gap near the Dirac points leading 
to the quantum- anomalous Hall effect. We also find that the noble metal atoms (i.e. Cu, Ag, Au) 
prefer the top adsorption site, and the dominant inequality of the AB sublattice potential opens 
another kind of nontrivial bulk gap exhibiting the quantum- valley Hall effect. 

PACS numbers: 73.22.Pr, 75.50.Pp, 75.70.Tj 



Introduction — . After an initial rush [1] since its 
first experimental exfoliation [2], graphene research has 
turned toward the application of graphene-based elec- 
tronics [l|,[3[. Graphene itself is a zero-gap semiconductor 
characterized by the linear Dirac-type dispersion, which 
is closely related to many amazing properties, e.g. high 
electron mobility and half-integer quantum Hall effect 
(per spin and valley). However, a band gap in graphene 
is highly desirable for designing semiconductor devices. 
So far, there are several proposals of engineering a band 
gap, such as by employing a staggered AB sublattice po- 
tential [4 , strain effect [5j] , intrinsic spin-orbit coupling [gl 
in single layer graphene, or applying a perpendicular elec- 
tric field in bilayer graphene [7|. 

Graphene is also considered as a promising candi- 
date for spintronics, which generally require imbalanced 
spin up and down carrier populations. However, the 
pristine graphene is nonmagnetic, and one need to em- 
ploy external methods to magnetize it, e.g. by decorat- 
ing hydrogen, sulfur or metal adatoms on graphene [H- 
[lO'], or substituting metal impurities for carbon atoms 
in graphene [11]. Alternatively, spin-orbit coupling is 
thought to be another source for manipulating the spin 
degrees of freedom through electric means [12]. In 
graphene there are two kinds of spin-orbit coupHngs: in- 
trinsic and extrinsic [6[. The former one has been shown 
to be unrealistically weak [is]. The latter one, which 
is also known as Rashba spin-orbit coupling, arises from 
the top-bottom layer symmetry breaking due to the pres- 
ence of a substrate or strong perpendicular electric field, 
which was reported to be remarkably large in graphene 
placed on top of Ni(lll) surface [14]. 

In this Letter, we present a systematic investigation 
on the magnetic and electronic properties of graphene 
adsorbed with Sd-transition and noble metal atoms in- 
cluding the spin-orbit coupling using first principles cal- 
culation methods. The results show that most transi- 




FIG. 1: (Color online) (a) A 4 x 4 supercell of graphene. ai,2 
(a| 2) indicates the primitive vectors for the 1x1 (4x4) su- 
percell. Three possible adsorption sites of single adatom on 
graphene are labeled as: hollow (H), top (T), and bridge (B). 
(b) Corresponding reciprocal momentum space structures: 
red and black vectors are for the 1x1 and 4x4 supercells, 
respectively. 



tion metal adatoms favor the hollow site, while the noble 
metal adatoms prefer the top site. We find that both 
types of adsorption give rise to nontrivial bulk gaps. Fur- 
ther Berry-phase analysis [ig] demonstrates that the gap 
opening mechanisms are completely distinct, i.e. one 
arises from the joint effect of exchange field and Rashba 
spin-orbit coupling leading to the quantum- anomalous 
Hall (QAH) effect, while the other one originates from 
the inequality of the AB sublattice potential resulting 
in the quantum- valley Hall (QVH) effect. Our findings 
should not only provide a new scheme for the gap open- 
ing in graphene, but also show great potential for the 
realization of the long-sought quantum anomalous Hall 
effect and the dissipationless valley tronics. 

Computational Methods — . The adatom-graphene 
system is modeled using one metal atom on top of a 
N X N (A^=3, 4) supercell of graphene. As illustrated in 
Fig. [11(a), a supercell is composed of 2N'^ carbon atoms 
and one metal atom, and the single atom has three possi- 
ble adsorption sites: hollow (H), top (T), and bridge (B). 
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TABLE I: Energetic and structural properties of 12 metal atoms located at the lowest energy adsorption sites on the 3x3 
and 4x4 supercells of graphene. The properties listed include the favored adsorption sites, adatom-graphene distance d (A), 
adsorption energy 6E (in units of eV/aiom), bulk gap A (meV), charge transfer SQ (e) from adatoms to graphene, and the 
magnetization of the adatoms /Jag (i^b) on the 4x4 supercell. Note that data inside the parentheses are used for comparison 
at the hollow site. 0.0* means that the global band gap is zero, but the local gap is finite around K/K' points for the 4x4 
supercell or V points for the 3x3 supercell. 



3x3 supercell 4x4 supercell 



Atoms 


site 


d 


^E 


A 


5Q 


site 


d 


^E 


A 


8Q 




Sc 


H 


L92 


1.26 


2.5 


0.97 


H 


1.97 


1.43 


2.5 


0.89 


2.30 


Ti 


H 


L86 


1.76 


0.0* 


0.88 


H 


1.86 


1.75 


1.4 


0.87 


3.41 


V 


H 


L89 


0.94 


0.0* 


0.67 


H 


1.87 


1.01 


2.2 


0.71 


4.51 


Cr 


B(H) 


2.33(2.17) 


0.14(0.13) 


0.0(0.0*) 


0.35(0.41) 


B(H) 


2.37(2.15) 


0.16(0.15) 


0.0* (2.5) 


0.37(0.45) 


5.64(5.52) 


Mn 


H 


2.07 


0.14 


4.5 


0.56 


H 


2.09 


0.14 


4.5 


0.57 


5.42 


Fe 


H 


1.56 


0.98 


0.0* 


0.81 


H 


1.56 


0.95 


5.5 


0.81 


1.97 


Co 


H 


1.55 


1.31 


0.0* 


0.63 


H 


1.53 


1.67 


0.0* 


0.66 


0.99 


Ni 


H 


1.58 


1.37 


60 


0.51 


H 


1.58 


1.39 


0.0 


0.52 


0.00 


Cu 


T(H) 


2.24(2.17) 


0.24(0.10) 


0.7(0.4) 


0.16(0.14) 


T(H) 


2.25(2.17) 


0.25(0.10) 


38.8(7.0) 


0.15(0.15) 


0.90(0.92) 


Zn 


H 


3.78 


0.03 


0.6 


0.02 


H 


3.78 


0.02 


0.0 


0.02 


0.00 


(LDA) 


H 


3.04 


0.14 


13.9 


0.03 


H 


3.02 


0.14 


0.0 


0.03 


0.00 


Au 


T(H) 


3.13(3.47) 


0.14(0.13) 


0.0* (2.9) 


-0.08(-0.09) 


T(H) 


3.45(3.54) 


0.16(0.16) 


0.6(1.7) 


-0.02(-0.15) 


0.84(0.84) 


(LDA) 


T(H) 


2.37(2.43) 


1.20(0.94) 


9.4(2.0) 


-0.06(-0.05) 


T(H) 


2.40(2.52) 


0.77(0.50) 


161.4(8.6) 


-0.09(-0.09) 


1.03(0.96) 


Ag 


T(H) 


3.57(3.69) 


0.03(0.03) 


0.0(1.5) 


0.04(0.02) 


T(H) 


3.49(3.58) 


0.03(0.03) 


3.0(2.3) 


0.02(0.02) 


1.00(1.00) 


(LDA) 


T(H) 


2.45(2.52) 


0.82(0.76) 


2.8(6.8) 


0.08(0.07) 


T(H) 


2.44(2.51) 


0.38(0.31) 


21.8(8.7) 


0.11(0.08) 


0.87(0.93) 



Panel (b) plots the Brillouin-zone of the 4x4 supercell 
comparing to that of the 1x1 supercell with well-defined 
valley indices K and K' . 

The first principles calculations were performed us- 
ing the projected- augmented- wave method as im- 
plemented in the Vienna Ab-initio Simulation Pack- 
age (VASP) [18]. The generalized gradient approxi- 
mation (GGA) exchange-correlation functional [19j was 
mainly used except those specified by the local density 
approximation (LDA) in Table HI The kinetic energy cut- 
off was set to be 500 eV, and the experimental lattice con- 
stant of graphene a = 2.46 A was used. During the struc- 
ture relaxation, all atoms were allowed to relax along the 
normal direction of graphene and all parameters were 
chosen to converge the forces to less than 0.01 eV/A. 
The first Brillouin-zone integration was carried out by 
using the 12 x 12 x 1 and 6x6x1 Monkhorst-Pack grids 
for the 3x3 and 4x4 supercells, respectively. A vacuum 
buffer space of 15 A was set to prevent the interaction 
between adjacent slabs. 

Adsorption Analysis — . Table |T] summarizes the sta- 
ble adsorption sites, energetic and structural properties 
of 12 kind of 3(i-transition and noble metal atoms for 
the 3x3 and 4x4 supercells. The stable adsorption 
site, distance d and magnetization of adatoms jiag were 
obtained by only considering the magnetization, while 
other quantities were calculated by further including the 
spin-orbit coupling. The adsorption energy is defined 
as: 6E = Ea Eg — Eag, where Ea, Eg and Eag are 



the energies of the isolated atom, N x N supercell of 
graphene, and adatom-A/" x N supercell of graphene, re- 
spectively. The charge transfer 5Q was calculated based 
on the Bader charge analysis [20]. For clear demonstra- 
tion, we divide the adsorption analysis into two parts: 

(a) For the 3(i-transition metal adatoms, we observe 
that most of them except Cr favor the hollow site, (i) 
For Sc, Ti, V, Fe, Co, and Ni, the adsorption distances 
are very short (less than 2 A), and the adsorption energy 
SE and the charge transfer are very large, which together 
indicate a strong hybridization; (ii) For Cr and Mn, the 
adsorption distance is slightly larger than 2 A, the ad- 
sorption energy is relatively weak but the charge transfer 
is still remarkably large, which corresponds to a moder- 
ately weak adsorption; (iii) For Zn, both the large ad- 
sorption distance and extremely small adsorption energy 
and charge transfer signal an extremely weak hybridiza- 
tion. In the following, we show that the reconstruction of 
the electronic configurations of the adsorbed metal atoms 
is intimately dependent on the strength of adsorption on 
the graphene sheet. 

In Table HI we notice that the magnetic moments of 
adatoms jUag are greatly altered compared to that of the 
isolated atoms. This is found to obey the following rules: 
(1) When the Sd-shell is less than half filled (i.e. Sc, 
Ti, V), the strong hybridization lowers the energy of the 
3(i-orbital. This makes the 4s electrons transfer to the 
unoccupied 3(i-orbital. Since the resulting 3(i-shell is not 
over half-filled, the Hund's rule dictates that all the 3d 
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electrons possess the same spin-polarization leading to 
the increase of the magnetic moments /j^ag] (2) When 
the 3(i-shell is exactly half-filled (i.e. Cr, Mn), the ad- 
sorption is weak therefore the 3(i-orbital energy is less 
affected during the hybridization. Based on the Hund's 
rule, all 3(i-electrons are equally spin-polarized showing 
a maximum magnetization. In the infiuence of the high 
magnetization from the 3(i-shell, the 4s electrons tend to 
align with the spin of the 3d electrons, which slightly in- 
creases the magnetic moment; (3) When the 3(i-shell is 
over half-filled (i.e. Fe, Co, Ni), the strong hybridization 
lowers the 3(i-shell energy and the transfer of 4s electrons 
to 3(i-orbital decreases the magnetic moment due to the 
occupying of the unpaired 3(i-orbital [9]. In particular, 
for Ni the two 4s electrons are transferred to 3(i-orbital 
and form a closed 3d shell, giving rise to a vanishing mag- 
netic moment; (4) For atoms with closed 3d and 4s shells 
(i.e. Zn), it behaves like a inert atom leaving graphene 
nearly unaffected. 

(b) For the noble metal adatoms, we find that they are 
stable at the top adsorption site. Consistent with Ref. [sl, 
the long adsorption distance, low adsorption energy, and 
weak charge transfer originate from the physical adsorp- 
tion involving the van de Waals forces. In Table HI we 
also provide the LDA results for Au and Ag adsorption. 
The major differences are that the adsorption distance d 
is dramatically decreased and the adsorption energy 5E 
is exponentially increased. Because of the closed d shells, 
the single s valence electron contributes to the magnetic 
moment jj^ag ^ IfiB- 

Due to the magnetic proximity effect with the magnetic 
adatoms, graphene can be magnetized. In the absence of 
spin-orbit coupling, we estimated the induced exchange 
field [Til near the Dirac points to be A=145 (Sc), 277 (Ti), 
199 (V), 359 (Cr), 181 (Mn), 57 (Fe), 462 (Co), (Ni), 
51 (Cu) in units of meV for the hollow-site adsorption on 
the 4x4 supercell. 

Band Structures and Berry Curvatures — . We classify 
the twelve atoms into three groups: (1) atoms (e.g. Sc, 
Mn, Fe) adsorbed at the hollow site, and graphene is 
magnetized; (2) atoms (e.g. Ni) adsorbed at the hollow 
site, but graphene is not magnetized; (3) atoms (e.g. Cu) 
adsorbed at the top site. To capture the physics arising 
from the Dirac points, we center on the adsorption on the 
4x4 supercells. The adsorption on the 3x3 supercells 
will be discussed in the supplementary materials. 

Figure [2] plots the band structures of adatom-graphene 
for (a) Sc, (b) Mn, (c) Fe, (d) Cu, (e) Ni adsorbed at 
the hollow sites and (f) Cu adsorbed at the top site. 
Panels (ai)-(fi) exhibit the full bulk bands along the high 
symmetry lines, and panels (a2)-(f2) zoom in the bands 
near the Dirac points K and K' . We see that sizable 
bulk band gaps open at the Dirac points in panels (a2)- 
(d2) and (f2). The presence of such a bulk gap indicates 
an insulating state. However, we note that all the Fermi 
levels lie outside the bulk gaps except that in panel (C2) 




FIG. 2: (Color online) Panels (ai)-(ei): Full bulk band struc- 
tures based on GGA for Sc (ai), Mn (bi), Fe (ci), Cu (di), 
Ni (ei) adsorbed at the hollow site on a 4 x 4 supercell of 
graphene along high symmetry lines by including the spin- 
orbit coupling. Panel (fi): Full bulk band structure based 
on GGA for Cu adsorbed at top site on a 4 x 4 supercell of 
graphene. Panels (a2)-(f2): Zooming in of the bands around 
K and Dirac points (in black), and the Berry curvatures 
Q{k) for the whole valance bands (in red). 



for Fe-adsorption. Thus, to realize this insulator one has 
to artificially adjust (e.g. by applying a gate voltage) the 
Fermi levels to be inside the gap. 

To further identify the nontrivial topological properties 
of the obtained insulators, we turn to the Berry phase 
analysis in the momentum space [l6|. By constructing 
the Bloch wavefunctions tjj from self-consistent poten- 
tials, the Berry curvature fl{k) can be obtained using 
the following formula : 

n 

n (U\ - 2Im ( V^nfc \Vx\i^n'k){'^n'k\Vy\ j^nk ) x 
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where the summation is over ah n bands below the bulk 
band gap, fn is the Fermi-Dirac distribution function, 
ujn = En/Ti^ and v^f^y^ is the velocity operator. In pan- 
els (a2)-(d2) and (£2), we also plot the Berry curvatures 
Vt{k) in dotted curves. One can observe that in panels 
(a2)-(d2) the Berry curvatures Vt K point share the 
same sign as that at K' point. By integrating the Berry 
curvatures over the first Brillouin zone, it gives rise to a 
nonzero Chern number C = +2 for the insulating state 
of Sc, Fe, Cu adsorption, and C = — 2 for the insulating 
state of the Mn adsorption. The non-vanishing Chern 
number C in the absence of magnetic field signifies a QAH 
state [24], which was proposed theoretically [5] but still 
not realized in experiment. In the subsequent section, we 
show the mechanism of the formation of the QAH effect 
in graphene. 

Figure [2] clearly shows that the adatoms of group (i) 
can open a non-trivial bulk band gap and achieve the 
QAH effect, but that of group (ii) can not. There- 
fore, we conclude that the magnetization in graphene, 
which breaks the time-reversal symmetry, is an essen- 
tial ingredient for the QAH effect. Alternatively, we find 
that no bulk band gap opens for the group (i) atom- 
adsorption when the spin-orbit coupling is switched off. 
This makes the spin-orbit coupling the other key ingre- 
dient of the QAH effect. From Table HI one can notice 
that there is substantial charge transfer from the metal 
adatoms to graphene. Due to the asymmetry of the 
adatom-graphene structure, a sizable potential gradient 
is formed at the interface, which induces a considerable 
Rashba spin-orbit coupling Q. Panel (d2) for Cu ad- 
sorbed at the hollow site further supports that the 4s 
spin-polarized electron transfer can also realize the QAH 
state. Therefore we summarize that the realization of 
the QAH state is only dependent on the spin-polarized 
charge transfer, but independent of the details in the ad- 
sorption (i.e. strong/weak hybridization). Although the 
induced Rashba spin-orbit coupling and magnetization 
in graphene are non-uniformly distributed in the struc- 
ture, the gap opening mechanism near the Dirac points is 
similar to that of the tight-binding model with uniformly 
distributed Rashba spin-orbit coupling and magnetiza- 
tion [2^: first, the magnetization lifts the spin up and 
down bands, accompanying with the bands crossing; sec- 
ond, the Rashba spin-orbit coupling opens a bulk gap at 
the bands crossing points. 

Because of the Van de Vaals force, the noble metal 
adatoms favor the top adsorption site. Besides the result- 
ing Rashba spin-orbit coupling and exchange field, the 
major difference from the hollow-site adsorption is the 
extra introduction of the inequality of AB-sublattice po- 
tential, which breaks the spacial inversion symmetry. For 
example, we have considered the case with Cu adsorbed 
at the top site as shown in the panel (f2). Through plot- 
ting the Berry curvature distribution, we observe that Vt 
at K point exhibits opposites sign as that at K' point. 



This strongly signals a QVH state with canceling Chern 
numbers [2J], i.e. C = 0, but Ck = —Ck' = 1. We can 
attribute this result to the leading role of the AB sublat- 
tice inequality suppressing the joint effect from Rashba 
spin-orbit coupling and exchange field. 

Conclusions — . We study the magnetic and electronic 
properties of graphene adsorbed with 12 metal atoms us- 
ing first principles calculation methods. We show that 
the QAH effect can be realized by doping the magnetic 
adatoms (i.e. Sc, Ti, V, Mn, Fe, Cu, Ag, and Au) at 
the hollow site of graphene, which originates from the 
joint effect between the Rashba spin-orbit coupling and 
magnetization. We also find that when the noble metal 
atoms (i.e. Cu, Ag, and Au) are adsorbed at the top 
site, the QVH effect can be observed due to the inversion 
symmetry breaking from the inequality of the AB sublat- 
tice potential. Our findings of the QAH and QVH states 
in the adatom-graphene systems not only provides a new 
gap opening scheme for the industrial application, but 
also pave the way for the realization of the dissipationless 
charge/ valley current in spintronics and valley tronics. 
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Supporting Materials 

In Section I, we shall discuss the band structures of the 
3x3 supercells of graphene adsorbed with various transi- 
tion metal atoms. Our results confirm that the quantum 
anomalous Hall (QAH) state can also be formed even 
though the non-trivial bulk gap opens at the F point 
but not at the K/K' points in the 3x3 supercells of 
graphene. In Section II, we introduce the definitions of 
the QAH effect and quantum valley Hall (QVH) effect. 

Band structures for adatom-3 x 3 supercells of 
graphene 

In the main text, we have shown the band structures 
for the 4x4 supercells including the spin-orbit coupling 
and magnetization using the first-principles calculation 
method. The corresponding K and K' valleys are well- 
separated and distinguishable. We have shown that the 
resulting topological (QAH/QVH) states can be formed 
due to the non-trivial bulk gap opening around the K 
and K' points. 

For the the adsorption in the 3x3 supercells of 
graphene, due to the folding of both K and K' into the 
same F point, valley indices are indistinguishable. There- 
fore, it has not attracted so much interest as those struc- 
tures (i.e. 4x4, 5x5, 7x7) with distinguishable valleys. 
However, in experiment it may form the 3x3 super- 
cell structures during the adsorption of diluted atoms on 
top of graphene. Therefore, it is necessary to know that 
whether the topological states (i.e. QAH state) proposed 
in the 4x4 supercells of graphene can still survive or not 
in the 3x3 supercells of graphene. 



In Fig. [3] (a) we schematically plot the 3x3 supercell 
of graphene compared to the 1x1 supercell of graphene, 
with ai^2 and a\ 2 representing the primitive vectors for 
the 1x1 and 3x3 supercells. Panel (b) illustrates the 
Brillouin zone of the 3x3 supercell comparing to that 
of the 1x1 supercell. One can note that the original 
valleys K and K' are folded into the same F point in the 
3x3 supercell of graphene, which makes the two valleys 
indistinguishable . 

Figure 2] exhibits the band structures of 3 x 3 super- 
cells of graphene adsorbed with (a) Sc, (b) Mn, (c) Fe, 
(d) Cu, (e) Ni and (f) Co adatoms on top of the hollow 
position without involving spin degrees of freedom, i.e. 
nonmagnetic calculation. In the first column [i.e. pan- 
els (ai)-(fi)], we plot the full band structures along the 
high symmetry lines. From the band structures shown 
in the second column [i.e. panels (a2)-(f2)], the zooming 
in of the circled bands in the first column, we observe 
that large band gaps can open at the F point due to 
the inter-valley scattering between K and K' valleys, i.e. 
A = 71.1,190.6, 100.8, 37.7, 106.7, 119.4 in units of 
meV for the 3x3 supercell of graphene adsorbed with 
Sc, Mn, Fe, Cu, Ni, and Co, respectively. This is com- 
pletely distinct from that of the adsorption on the 4x4 
supercells of graphene in the main text, where the per- 
fect Dirac-cone dispersion holds at both K and K' points. 
In our calculations we have used the Hydrogen atom to 
saturate the dangling bond of Mn atom, because in the 
Mn-adsorption the dangling d-band is properly located 
inside the resulting bulk gap. In the first column one can 
see that there are two bulk band gaps opened at the F 
points. As shown in Fig. [5] the projected band analysis 
points out that the two gaps are essentially equivalent by 
combining the same orbital components of P^, dxz and 
dyz (or P^, dj:y and dx2_y2). 

Similar to the adsorption on the 4x4 supercell of 
graphene, when the magnetization is included, the spin 
up (down) bands are upward (downward) shifted (see 
Figure [6]). This results in three possible results: (i) As 
plotted in panels (a2)-(c2) and (f2) for Sc, Mn, Fe, and 
Co adsorption, the induced magnetization is so large (i.e. 
A=0.213, 0.432, 1.080, 0.442 in units of eV) that the orig- 
inal band gaps from the inter- valley scattering are closed 
due to the crossover between the spin-down bands from 
the conduction bands and the spin-up bands from the 
valence bands, (ii) For panel (d2) of Cu- adsorption, the 
magnetization is moderately large (i.e. A=0.035 eV) that 
the band gaps are just slightly decreased but not closed, 
(iii) For the Ni- adsorption in panel (e2) the band struc- 
ture is exactly the same as that shown in Fig. |4] (e2) 
because there is no magnetization of graphene induced 
from the Ni adsorption. 

When the spin-orbit coupling is further considered, we 
find that non-trivial band gaps can open at the cross- 
ing points between the spin-up and spin-down bands [see 
panels (a2)-(c2) and (f2) in Fig. [7], which shows similar 
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FIG. 3: (Color online) (a) A 3 x 3 supercell of graphene. ai,2 
(ai 2) indicates the primitive vectors for the 1x1 (3x3) 
supercell. Three possible adsorption sites of single adatom on 
graphene are labeled as: hollow (H), top (T), and bridge (B). 
(b) Corresponding reciprocal momentum space structures: 
61,2 and 6f 2 are for reciprocal vectors for the the 1x1 and 
3x3 supercells, respectively. K and points for the 3x3 
supercell of graphene are folded into the F point. 



origin as that in the adsorption on the 4x4 supercells 
in the main text. However, we find that only the bulk 
band gaps in panels (ai-bi) are globally opened for Sc 
(A = 2.5 meV) and Mn (A = 4.5 meV), while that for 
Fe and Co in panels (ci) and (fi) are just locally opened. 
Therefore in the following, we will only focus on the non- 
trivial topological properties for Sc and Mn adsorption. 
Through calculating the Berry curvatures Q along the 
high symmetry lines for all the valence bands below the 
bulk gap, we find that all the Berry curvatures ft share 
the same sign as shown in panels (a2) and (b2), resulting 
in the non- vanishing Chern number C. By integrating all 
the Berry curvatures below the band gaps, our calcula- 
tion shows that C = +2/ — 2 for the Sc/Mn adsorption. 
This again signals a QAH insulating state for the ener- 
gies inside the bulk gap which will be explained in detail 
in Section IL Since the Fermi energy is not located inside 
the band gap, one has to apply external gate voltage to 
adjust the Fermi level inside the band gap to realize the 
insulating state. While for the other two adsorption of 
Cu and Ni shown in panels (d2)-(e2), we can see that the 
spin-orbit coupling further decreases the bulk gaps. 

In summary, using the first-principles calculation 
methods, we find that for the 3x3 supercell of graphene 
adsorbed with atoms (i.e. Sc, Mn, Fe, Cu, Ni) on the 
top of hollow position, trivial bulk band gaps can open 
at the r point due to the inter- valley scattering without 
considering the spin degrees of freedom. For some partic- 
ular adatoms (i.e. Sc, Mn) when the magnetization and 
spin-orbit coupling are involved, we find that non-trivial 
bulk band gaps can open around the F point to realize 
the QAH effect. This finding not only generalizes our 
proposal of QAH effect in the 4x4 supercell to the 3x3 
supercell of graphene, but also enhances the experimental 
realizability of the first QAH effect in graphene. 



Definitions of Quantum Anomalous Hall and 
Quantum valley Hall effects 

Quantum Anomalous Hall Effect 

The breaking of the time-reversal symmetry is a nec- 
essary condition for the formation of the Hall effects. 
In general, there are two kinds of methods breaking the 
time-reversal symmetry: the external magnetic field and 
the intrinsic magnetization of the ferromagnetic materi- 
als. When the two-dimensional electron gas is subjected 
to a strong perpendicular magnetic field, Laudau-levels 
can be formed to exhibit the integer quantum Hall effect 
with the Hall conductance exactly quantized in units of 
e'^/h [ij. The precise quantization of Hall conductance 
can be understood in terms of a topological invariant 
known as the Chern number C, which must be an inte- 
ger Q . Therefore the integer quantum Hall conductance 
can be expressed as: 

(Txy=C— (2) 

Contrary to the ordinary Hall effect due to the mag- 
netic field, the intrinsic magnetization induced Hall ef- 
fect is termed as anomalous Hall effect. Though this 
anomalous Hall phenomenon has been known for over 
one century, its origin is still unclear [sl, |^ . The quan- 
tized version of the anomalous Hall effect, also named as 
quantum anomalous Hall (QAH) effect, was first theoret- 
ically proposed by Haldane in a Honeycomb toy model 
after the experimental observation of the integer quan- 
tum Hall effect in 1980s However, the QAH effect is 
still not observed in experiment. 

In our work on the metal adatom-graphene structures, 
we find that the QAH state can be observed in both 3x3 
and 4x4 supercells of graphene. In the main text the bulk 
gaps are opened at both K and K' points when the mag- 
netization and spin-orbit coupling are considered. Be- 
cause the Berry curvatures near K and for the 4x4 
supercell show the same sign, therefore the summation 
of the Berry curvatures of the whole Brilloin-zone is fi- 
nite, which indicates that the Chern number must be a 
nonzero integer. Since we do not apply any external mag- 
netic field, therefore we can conclude this non-vanishing 
Chern number corresponds to the QAH state. Surpris- 
ingly, in the Section I of this supplementary materials we 
further find that the adsorption on the 3x3 supercell 
can also contribute to the QAH effect through opening a 
nontrivial bulk band gap near the F point by involving 
the magnetization and spin-orbit coupling. 

Quantum Valley-Hall Eff*ect 

Quantum valley Hall (QVH) effect is named after the 
integer quantized but canceling quantum-Hall responses 
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at K and K' valley points. This effect can be arisen from 
the inversion- symmetry breaking, i.e. applying staggered 
AB sublattice potentials on the single layer graphene [6| , 
or considering gated bilayer graphene [7|]. In terms of 
Chern number representation, Ck = —Ck' 7^ but the 
total Chern number of the system is C = Ck + Ck' = 
0. One has to note that there is no edge modes for the 
Q VH effect in the single layer graphene ^, ^] , while for 
the bilayer graphene there is indeed gapless edge states 
inside the bulk band gap [9] . These edge states are robust 
against smooth-type disorders, which is protected by the 
large momentum separation between two valleys. 

In our studied adatom-3 x 3/4 x 4 structures, no QVH 
state can exist in the 3x3 supercell systems because of 
the inter- valley scattering which makes valleys indistin- 
guishable. In the 4x4 supercell with noble atom ad- 
sorbed graphene systems, our calculation shows that the 
Berry curvatures at K and K' are non- vanishing and op- 
posite, which indicates that valleys K and K' correspond 
to two opposite nonzero Chern numbers. Based on the 
definition of QVH effect, we can claim that our finding of 



the noble adatom adsorption on top of the 4x4 supercell 
leads to the QVH effect. 
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FIG. 4: (Color online) First column (ai-fi): full band structures of adatom-3 x 3 supercell of graphene for Sc, Mn, Fe, Cu, Ni, 
Co adsorbed on top of the hollow position without including magnetization and spin-orbit coupling along the high symmetry 
line. Second column (a2-f2): Zooming in of the circled bands around F point in panels (ai-fi). Bulk band gaps are opened at 
the F point due to the inter- valley scattering effect. Note that there are two equivalent band gaps at F point except in the 
band structure of Cu-adsorption. 




FIG. 5: (Color online) Projected bands analysis for the bands near the two bulk band gaps in Figure |4jbi). We can see that 
the bands near the two bulk gaps possess the same orbital components. 
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FIG. 6: (Color online) First column (ai-fi): full band structures of adatom-3 x 3 supercell of graphene for Sc, Mn, Fe, Cu, 
Ni, Co adsorbed on top of the hollow position by including the magnetization along the high symmetry line. Second column 
(a2-f2): Zooming in of the circled bands around F point in panels (ai-fi). Colors are used to distinguish the spin- up (black) and 
spin-down (red) bands. In the second column all the spin-up and spin-down bands are crossing due to the large magnetization 
except that for Cu and Ni adsorption. In panel (d2) for Cu- adsorption the magnetization is less than the original band gap 
amplitude, therefore the relative shift between the spin-up and spin-down bands does not cross to close the band gap. In panel 
(e2) for Ni-adsorption the band structure is exactly the same as that in Figure [41(e) because of the vanishing magnetic moment 
of Ni-adatom. 
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FIG. 7: (Color online) First column (ai-fi): full band structures of adatom-3 x 3 supercell of graphene for Sc, Mn, Fe, Cu, 
Ni, Co adsorbed on top of the hollow position by including both the magnetization and spin-orbit coupling along the high 
symmetry line. Second column (a2-f2): Zooming in of the circled bands around F point in panels (ai-fi). In panels (a2-C2) and 
(f2), nontrivial bulk band gap open at the crossing points. However, only bulk band gaps are globally opened in panels (a2-b2) 
for Sc, Mn adsorption, while the gaps in panels (C2) and (f2) are just locally opened. Therefore only the Sc/Mn adsorbed 3x3 
graphene structure can exhibit the nontrivial insulating properties when the Fermi-level is artificially adjusted to the bulk band 
gaps. In addition, we also plot the total Berry curvatures Q(k) in panels (a2-b2) for all the valence bands below the bulk band 
gaps. 



